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Abstract. We study the primordial non-Gaussinity predicted from simple models of inflation 
with a linear potential and superimposed oscillations. This generic form of the potential is 
predicted by the axion monodromy inflation model, that has recently been proposed as a 
possible realization of chaotic inflation in string theory, where the monodromy from wrapped 
branes extends the range of the closed string axions to beyond the Planck scale. The 
superimposed oscillations in the potential can lead to new signatures in the CMB spectrum 
and bispectrum. In particular the bispectrum will have a new distinct shape. We calculate 
the power spectrum and bispectrum of curvature perturbations in the model, as well as make 
analytic estimates in various limiting cases. From the numerical analysis we find that for a wide 
range of allowed parameters the model produces a feature in the bispectrum with /nl ^5 — 50 
or larger while the power spectrum is almost featureless. This model is therefore an example of 
a string-inspired inflationary model which is testable mainly through its non-Gaussian features. 
Finally we provide a simple analytic fitting formula for the bispectrum which is accurate to 
approximately 5 % in all cases, and easily implementable in codes designed to provide non- 
Gaussian templates for CMB analyses. 
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1. Introduction 

In recent years there has been an increasing interest in the primordial non-Gaussianity predicted 
from inflation. Single field slow-roll inflation generally predicts non-Gaussianities to be 
unobservably small with present day techniques [IH3]. However, in non-minimal models of 
inflation the non-Gaussian signatures can be observable and non-Gaussianity is therefore an 
interesting observational window into the complexity of inflation. 

Both from an observational and theoretical point of view the effective model of inflation is 
still uncertain. Since we do not know the microscopic origin of the effective theory of inflation, 
then the simplest theory of single field slow-roll inflation seems preferred from an Occam's razor 
point of view. Thus, with no prejudice about the underlying microscopic theory of inflation, we 
do not expect large non-Gaussianities. On the other hand, when inflation is embedded into a 
fundamental theory like string theory, it typically carries with it more complicated structures, 
and therefore from a model building point of view more complicated models of inflation can be 
attractive. 

The simplest possible model of inflation, which is a single field model of chaotic inflation 
with a monomial potential, requires the field values of the inflaton to be trans-Planckian. This 
means that higher dimensional operators has to be prohibited by some symmetry in order not 
to destroy the flatness of the potential. In Natural Inflation models OE] this is achieved by 
introducing an approximate shift symmetry, and assuming that the inflaton is the pseudo- 
Nambu-Goldstone boson mode (axion) of the theory. It has generally been considered non- 
straightforward to realize this idea in string theory, because the field values would typically 
be restricted to sub-Plackian values. This has led to the proposal that one could extend the 
field range by having many axions [7J. Recently another proposal has appeared^, where the 
field range is extended by a so called monodromy mechanism [Hl[l2]- While not generic in 
this model, the model can have interesting oscillatory features embedded on top of the linear 
inflaton potential, depending on the details of the compactification [13J. 

Oscillatory features in the inflaton potential has already been studied by Chen, Easther 
and Lim as a way of generating large non-Gaussianity by a kind of resonance effect [H]. The 
axion monodromy inflation model is a concrete model, which could lead to an observational 
signature of this particular type of non-Gaussinity. Since the small oscillations in the potential 
will be very hard to detect directly in the power spectrum, the signature in the bispectrum 
may provide the best way of experimentally verify such a scenario. 

Non-Gaussianity typically manifests itself through a non- vanishing three-point correlation 
function of the curvature perturbation, and is generally quantised in terms of /jvl, which 
measures the deviation in the curvature perturbation ( from a Gaussian distributed variable ( g 

f Alternative proposals have also been promoted in [8HI0] 
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through [2J 

C = ( 9 - \InlQ ■ (i) 

At present, there are only good observational constraints on a scale independent f^L 
(f l N C L al ), which gives -10 < fftf < 70 at 95% C.L. [IMS] (note the preference for > 
which in some analyses is significant at more than 95% C.L.), but in the present model Jnl 
will be scale dependent, and it is not clear how the bounds on f 1 ^ 1 translate into a bound on 
the Jnl in the present model, but one presumingly loose some sensitivity by only constraining 
I'nT 1 experimentally. 

In the axion monodromy model, there are also non-trivial features in the power spectrum, 
and in the usual definition of /nl the power spectrum can therefore introduce some spurious 
scale dependence, which is not coming from the three-point function. Though this effect is at 
less than the per cent level for the specific models considered in this paper, we discuss a suitable 
alternative. 

In order to calculate quantitatively the bispectrum in the axion monodromy inflation 
model, we have developed a code for computing non-Gaussianity in general models of single 
field inflation with a minimal kinetic term. 



2. Bispectra in single field inflation 



In the case of single field inflation with some generic potential, one can write the minimally 
coupled Einstein-Hilbert action as 



S 



l -R-\{d<pY-v{<p) 



(2) 



At zeroth order, this leads to the ordinary background equation of motion for the scalar field 

(3) 



a a d(p 



and the time evolution of the slow-roll background can typically be expanded in the small 
slow-roll parameters 

# 1 2 
e = ~lP ~ 2~lP 

€ 



V = 



(4) 
(5) 



eH Hep 
where H = a/a is the Hubble parameter. 

Since we want to calculate correlation functions of fluctuations around this background, it 
is convenient to apply the ADM formalism [2D]. Here the perturbed metric takes the form 

ds 2 = -N 2 dt 2 + h i:j (dx i + N l dt) (dx j + N j dt) (6) 
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where the lapse and the shift N, N l acts as Langrange multipliers when the perturbed metric 
is inserted into the action. Only (f> and hij are dynamical variables, and in a convenient gauge 
that fix time and spatial reparametrizations, one can write 

hi, = a 2 (t) [e 2 % + 7*;] , (7) 

while having no fluctuations in the inflaton field. The physical degrees of freedom in this gauge 
are the scalar curvature perturbation £ and the transverse and traceless tensor mode 7^. We can 
safely ignore tensor perturbation in this work, as any contribution from tensor perturbations 
to the scalar three-function only arise through loop diagrams. 

Inserting the perturbed gauge fixed metric into the action, we can derive the action for 
the curvature perturbation iteratively to any order in perturbation theory. From the quadratic 
action of ( 

S= - [ dtd 3 x- 



,, /; L «V- «(<K)\ , (8) 

upon quantizing the ( field in the Bunch-Davies vacuum, one obtains the power spectrum, 
P^(k), of curvature perturbations 

(C kl Ck 2 ) = (27r) 3 5 3 (k 1 + k 2 )P c (fc) , (9) 

which is related to the usual scale independent power spectrum by P{(k) = (27r 2 ) / k 3 V^(k) . 
Similarly one can define the bispectrum Bq by 

(C(k 1 )C(k 2 )C(k 3 )) = (27r) 3 5(^k a ) J B c (k 1 ,k 2 ,k 3 ) , (10) 

a 

The strength of the non-Gaussian signal in the bispectrum can conveniently parameterised in 
terms of the scale independent non-linearity parameter /nl, by defining 
6 

B ( = -fNL[Pd k i) p <( k 2) + 2 permutations] . (11) 



Now, in order to calculate the three-point function one calculates the action for the 
curvature perturbations, £, to third ordeij], which, up to total derivatives, leads to the 
interaction Hamiltonian [2] 

H int (r) = ~ J d*x [ae 2 CC 2 + ae^dQ 2 - 2eC(d()(d X ) 

+?o/C 2 C + f (dC)(d x )(d 2 x) + Ud 2 Q{dx?] , (12) 

2 2a 4a J 

where d 2 x = eC 

One can, through this result, evaluate the three-point function of curvature perturbations 
after horizon exit, by means of the in-in formalism: 



(Cfc 1 (r)Cfa(r)Cfc 8 (r)> = -i / dra ([C kl (r)C k2 (r)C k3 (r), H mt (r')]) , (13) 

J TQ 

I The fourth order action has been calculated in [21], [3], [22] . 
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together with a field redefinition accounting for boundary terms that was neglected in eq. (TT21 

(Cfa(r)Cta(r)C*B(r)> = (Cfe(r)Cb(r)Cta(r))+^(S(r)&(r)i(r))+^ .(14) 

In a scenario where there are periodic ripples in the inflaton potential, the slow-roll 
parameters will generally oscillate. In such a scenario, it was demonstrated by Chen, Easther 
and Lim [14j that when the perturbations modes are inside the horizon and oscillate, they can 
interfere constructively with the oscillations in the couplings in the interaction Hamiltonian 
above, which are determined by the slow-roll parameters e, rf. As a consequence there can be 
a resonant amplification of non-Gaussianity from the bispectrum. 



3. The effective model of inflation from axion monodromy 

It has been demonstrated that in the large field limit, which is valid during inflation, the axion 
monodromy model is described by a linear potential for the axion <p and an axion action of the 
form |12] 



3$ = J d^x^fg ^-(<90) 2 — /i 3 0^ + corrections 



(15) 



One can show that in order to obtain 60 e-folds of inflation with the right level of curvature 
perturbations in the leading order linear potential, one should have 

0~llM pl (16) 
/i~6x l(T 4 M pl (17) 

as initial conditions. Though the brane-induced inflaton potential described above is the leading 
effect for breaking shift symmetries in the class of models considered here, there are other effect 
present as well. Among the more important is production of instantons, which, in this case, 
give rise to periodic corrections to the potential. One can therefore write down the effective 
potential containing these two contributions as [32] 



V eff = /rV(l + ai cos(0//)) + a 2 Mt cos(0//) , (18) 

where M s is the string scale and a, are dimensionless parameters which are expected to be much 
smaller than one, ctj -C 1, depending on the details of the string theory compactifications. 

We will study the two cases a.\ <C a 2 and cti ^> a 2 , leading to two limiting potentials. For 

Oil 3> « 2 

V#} = ^V(l + acos(0//)), (19) 

where we will assume a = a\ ~ 1CT 5 , such that effect on the power spectrum is at the per cent 
level and not yet excluded by data. The other case is a,\ a 2 where 

V%) =/x 3 + A 4 cos(0//), (20) 
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and where we will assume A 4 = a 2 M^ is such that the effect on the power spectrum is at the 
per cent level. 

The scalar power spectrum for the potential in the second model was computed analytically 
in [13] . We will return to this point later in Eq. ( I36lh 

3.1. Isosceles limit estimates of the bispectrum 

As discussed in section two, one can write f^L as a combination of the three-point function of 
the Fourier transform of ( and the power spectra of the curvature perturbations, yielding 

, = _io (hhh) 3 (CfciCfc 2 Cfc 3 ) (21) 

fNL 3 V{k x )V{k 2 )kl + 2 perms. (2vr) 7 5( 3 ) (fe a + k 2 + k 3 ) ' 1 ' 

where k% is the magnitude of the momentum fej. For the case of our models the fluctuations in 
the power spectrum is of the order one per cent. However if one should concern oneself with 
models with larger fluctuations in the power spectrum, one may benefit from using a similar 
definition, given by 

f = 10 (hk 2 k 3 f (Ck-XM , v 

JNL 3 V\k\ + k\ + k 3 3 ) (2tt)^(3) (k, + k 2 + k 3 )' { ] 

where V is approximately the amplitude of the power spectrum at the fc's considered. It should 
be noted that this definition reduces to the k ^ kz presented in [H] in the equilateral limit, but 
differs in the squeezed limit, where grows as ^ for triangles of the type k\ = k 2 k 3 . 

We proceed to find an estimate for the case k\ = k 2 > k 3 . Though such an estimate 
does not exist in general, we can obtain one by interpolating between the equilateral limit 
(fci = k 2 = k 3 ) and the squeezed limit {k\ = k 2 3> k 3 ). 

3.1.1. Equilateral Limit The model is essentially a linear slow- roll potential with a small 
oscillating perturbation. We can therefore follow the approach of Chen, Easther and Lim [TJ] 
and calculate the first order effect of the small oscillating term as a small correction. In this 
way one has e = eo + e osc , r\ = i] Q + r] osc . As the background solution (subscript 0), one can start 
from 

3#0 + fx 3 = (23) 



to obtain 



and 



2 2/3 



vV /2 t + 20, 3/2 ) 2 / 3 (24) 



1 2 

e o = 7T72 ' Vo = -jo ■ (25) 
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As argued in [H] one can estimate /tvl m the equilateral limit by 



(eq) 



In if 



+ phase 



(26) 



in the equilateral limit. Here K = k\ + k 2 + k 3 . The factor 2 missing, compared to [H], is 
due to the underlying linear potential here, as opposed to a quadratic in [14]. The resonance 
amplitude can be estimated as 



Ja 9 H, 

where tja is the amplitude of fj osc , which can be found from 



(27) 



(2i 



H osc A 4 1 

-7— e = -3 — — cos 

—^0 = -6^7 sin I — J 29 
A 4 -model 



We therefore find the resonance amplitude to be 
10 A 4 1 



(eq) 



9 0/i 3 /5/201/2 

A 4 -model 



(32) 



Hose ( 00 

— — e = —oa cos — 
H \f 



c osc 
Vosc = —VO 

eo 



c a0o . 
-o — y— sin 



./ 



a-model 



10 



-a- 



(eq) 

■1 ^ 0,^5/201/2 

a-model 



(30) 
(31) 



(33) 



3.1.2. Squeezed limit Let us now consider the different limit k\ = k 2 3> k% (for simplicity 
we define k\ = k 2 = k and k% = m). In this limit the long wavelength mode of £k 3 , will act 
as a constant rescaling of the background of the two shorter wavelength modes. As shown 
by Maldacena j2], one can thus calculate the three point correlation function in this limit by 
calculating the correlation of the long wavelength mode with the variation of the two-point 
function of the short wavelength modes on the background of the long wavelength mode 



lim ■ (CkiCkaCka) = (Ck 3 (CkiCk2>B) • ( 34 ) 
fc 3 ->0 

Then by a Taylor expansion of (CkiCk 2 ),B 011 the unperturbed background, and using d/d( —¥ 
kd/dk (since the effect of ( is to take k — >■ k — k(), one finds 



hm (C kl Ck 2 Ck 3 ) = -(27T) 3 5 3 (k 1+ k 2 + k3)Pc(fc3)^^ynW 



K - l)(2 7 r) 3 <5 3 (k 1 + k 2 + k 3 )P ( (k 3 )P c (k) (35) 
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which can be compared to eq. ffTTj) . and one can read off f$ = (5/12) (n s - 1). Here n s is the 
spectral index. 

As it was shown in [13] the scalar power spectrum for the potential in the second model 
can be written in the form 

h s — l 



V c (k) = V c (h) 



k 



1 + 5n s cos 



with the short hand notation of 1131: 



5n, 



12A 4 



- rnth 
8 LULU V 2/0, 



A 4 -model 



(37) 



./ 



8n, 



(36) 



12a</>* 



- rnth 
8 UJLU V 2/0, 



/ V 1 + (W 2 

a-model 



(3? 



The d>h is the field value, when the momentum k crosses the horizon. 



-21.il 1^1- <> - 



(39) 



The starred quantities are merely an arbitrarily chosen fix point, where k* is the momentum 
that crosses the horizon at = 0*. Though complicated, the expression for 5n s reduces to 



5n q w 6 



A 4 -model 



(40) 



5n, ~ 6 



2* 



/ 3/2 



a-model 



(41) 



for While the expression for d~n s is very precise for / > 10~ 3 it overestimates significantly 

the correct value for / <, 10~ 4 . 
As 5n s < 1 we find 



n a — 1 = rL 



1 + — — sm 



(42) 



to the first order in Sn s . With this one can find an expression for f^L i n the squeezed limit: 



f: 



(sq) 
NL 



/i S9) sin(^] . /, 



(»«) 



5 5n s 
12 74 



(43) 



to the first order in <5n s for k\ = k% ^> k 3 . 

This result can now be combined with eq. (126]) to form a general estimate for isosceles 
momentum triangles (k\ = k 2 > k^)\ 



SL~^(*.-l)+[fF + (f!f i -fF)j 



sm | + phase 



(44) 



where for completeness we have included the slow-roll contribution from the underlying linear 
potential, as the first term. 
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3.2. Semiclassical estimates of the trispectrum 

If the bispectrum can become large in a given model, one might be interested in the size of 
higher order correlation functions as well. While it is beyond the scope of this paper to study 
the trispectrum in details, we can still make some generic semiclassical estimates in certain 
limits. Most obviously we can study the trispectrum in the squeezed limit [21|25], but it is also 
possible to estimate a particular contribution to the non-Gassianity in the counter-collinear 
limit [3]. 

3.2.1. Squeezed limit One can imagine to make a definition of t/vl analogous to the definition 
of $nl in eq. ( !2T1) . relating it to the four-point function as 



l (hhhh) 3 (C(fci)C(fc 2 )C(fc 3 )C(fc 4 )) , . 

TNL ~v\2*Y £**? £ (3) (Ei kj) [b) 

By repeating the above calculation for k = ki = k 2 = k 3 ^> A; 4 under the same assumptions, 
one can get an estimate of the behaviour of the trispectrum in the squeezed limit 

lim (C kl Ck 2 Ck 3 Ck 4 ) = -(27r) 3 5 3 (k 1 + k 2 + k3 + k 4 )P c (A; 4 )l^^ (k 6 B c (k)) .(46) 

In order to evaluate this expression approximately, we can insert the estimate for the 
equilateral bispectrum obtained in eq. ( 126]) to find 



3 t{eq) V{kfV{k,) 

~ ^NL p 3 



tan ( ^ + $ 
2(n s - 1) + 



(47) 



However, we can obtain an even simpler estimate of the size of the trispectrum in specific 
configurations, by considering the double squeezed limit, where k = k\ = k 2 » k 3 » fc 4 . 
Iteratively, we then obtain 

lim (C kl Ck 2 Ck 3 Ck 4 ) ~ (n. - l) 2 (2 7 r) 3 5 3 (k 1 + k 2 + k 3 + k 4 )P c (A; 3 )P c (A ;4 )P c (fc048) 

k3,k4— >0 



This implies that in this limit tnl ~ (/] 



NL ) ■ 



3.2.2. Counter-collinear limit Let us consider the contribution to the four-point function 
from the exchange of a scalar mode between two pairs of external scalar modes. In the counter 
collinear limit where the momentum of the exchanged mode goes to zero k 4 + k 2 — > 0, the 
contribution to the four-point function from the exchange process can be expressed as the 
correlation of a pair of two-point functions (CkiCk 2 )> (Ck 3 Ck 4 ) due to the presence of the long 
wavelength scalar mode [1] 

, lim (C kl Ck 2 Ck 3 Ck4> = ((Ck 1 Ck 2 ) B (Ck 3 Ck 4 ) B ) 

ki+k 2 — >0 

= (n s - l) 2 (2vr) 3 5 3 (k 1 + k 2 + k 3 + ^)P c {k n )P C {ki)Pdh) ,(49) 
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with k\2 = |ki + k 2 |. 

Again from this contribution clearly t nl ~ (Inl) 2 - Thus if ~ 50, the find that there 
are large contributions to the trispectrum of order t^l ~ 2500, which are also in an interesting 
range for Planck, that should be sensitive to a ttvl ~ 560 [23] . 

4. Numerical results 

To validate our analytical estimates, and understand in which part of the parameter space the 
estimates are valid, we have developed a code, that can evaluate numerically the bispectrum 
for any single field inflation model. In the code we integrate the background evolution and up 
to three mode functions at a time, to solve for arbitrary momentum triangles. To integrate 
the different contributions to the bispectrum, we split the integrals in two: To account for the 
highly oscillating contribution inside the horizon, we use the method of Chen et al [TJJ, and 
partially integrate using the BD-vacuum as an approximation for the mode functions, while to 
evaluate the integral from typically 4 to 6 e-folds inside the horizon for the lowest wavelength 
mode until all modes are far outside the horizon, we integrate Eq. ( fT3l) with the full numerical 
solution to the mode functions. 

To compare the estimates to the numerical solution in a relevant part of the parameter 
space, we have chosen parameters that neither are ruled out by observations, nor by theoretical 
considerations, but still generates an interesting signal. As can be seen in Eqs. ( 126]) and (136]) . 
this model has the remarkable property that even though it is a single field model of inflation, 
it can generate a feature that is readily seen in the bispectrum, but not even with Planck, can 
be detected in the power spectrum. 

In Flauger et al p3] different observational and microphysical bounds of the model 
have been explored. They find that the non-detection of wiggles in the WMAP data limits 
A 4 //i 3 w acf) < 10 -4 , and that it is difficult to create realistic microphysical models that 
have / < 10 -4 . We have therefore chosen a range of models, specified in table H] that are 
characterised by having a currently non-observable wiggle amplitude power spectrum, but an 
/nl > 5. For the background potential we use \i ~ 6 x 10~ 4 M p i and the initial value of 
~ UMpi found in [12], that gives the correct overall amplitude in the power spectrum, and 
corresponds to approximately 60 e-folds before the end of inflation in the model. 

The power spectra of some of the models is shown in figure [1] and it can be seen that for 
the given values the oscillations in are at the per cent to per mill level. For some of the 
models the wiggles will not even be detected by the Planck satellite [26], while at the same 
time (see figure [I]) they have a significant amount of non-gaussianity. Ref. [26] studied the 
detectability of ripples due to trans-Planckian effects with Planck or a future cosmic variance 
limited experiment. In order for ripples to be unambiguously detected their frequency must 
lie within a range such that there is a number of oscillations in the observable fc-range and 
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Figure 2. The shape function x\ x\F{x\, x-x) as defined in Eq. ([50)) for the general ctg model. 
Note that by definition F is normalised to one for the equilateral form x± = x% = x$ = 1, 
and that we only plot unique triangles, i.e. triangles with x\ > X2 have been suppressed in the 
figure for clarity. 
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Figure 3. The bispectrum for isosceles triangles for the as and A§ models as a function of 
side ratio k/m together with the analytic model 



100.0 



10.0 



1.0 



0.1 




a data o 

A 4 estimate 

a estimate 



10 



100 



A 4 /(^ 3 f 5 ' 2 ),a/f 5/2 



1000 



Figure 4. Comparison of the analytic estimates of the amplitudes Ja of /nl to the numerical 
results. 
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4 x 10- 


-4 




5.95 x 10- 


-2 




general model 


0.1 


1 


3 x 10- 


-3 


1 x 10- 


-5 







5.63 x 10- 


-2 



Table 1. Numerically evaluated models. All models have fi ~ 6 x 10 M p i, and the comoving 
wave numbers A: are measured in units of (aH)o = (ai?)0=iiM p i- 



their frequency must be lower than the effective width of the experimental window function. 
With the frequency, u, defined such that P{k) oc sin(o;lnA;) [26] estimates that 1 <, w <, 50 
for it to be detectable. Even if it falls within this range the amplitude must of course also 
be sufficiently high. At approximately la the amplitude, A, should be larger than 0.0024 for 
Planck. Translating to our case we can make the identification u ~ 0*/ and A ~ 5n s . To take 
an example with the a-models: «3, «4, a§ have too small amplitude to be detected in P{k). a\ 
and «2 may be detectable, but a§ has u> ~ 50 and may have too fast oscillations for a detection. 
In conclusion, a number of the models with a measurably large non-Gaussianity would not have 
measurable wiggles in the power spectrum. 

We have used our code to probe a rather large range of parameter space, and in figure H] we 
show how the analytic estimates for the power spectrum and the bispectrum in the equilateral 
limit are in excellent agreement with the numerical results. Furthermore, we have also 
numerically evaluated squeezed triangles (see figure |3]), and found a simple functional form 
Eq. (jUJ) for general isosceles triangles, which agrees at the 5%-level with the numerical result. 

Given that the bispectrum in this model is nearly scale invariant, it makes sense to 
characterise the signal in terms of the shape function [27] 

k k 

F i x i = j-,%2 = Al) = fNL(k 1 ,k 2 ,k 3 )/f NL (k 3 ,k 3 ,k 3 ) (50) 
k 3 k 3 
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which measures the relative importance of different geometrical configurations compared to 
the equilateral configuration. In figure [2] it can be seen how the squeezed configurations are 
suppressed compared to the equilateral configuration. This is also observed in higher derivative 
models, and is in contrast to local type non-gaussianity from e.g. slow-roll inflation, where the 
squeezed limit is dominating the shape function [27]. But compared to other models in the 
literature, this models has the distinct feature of oscillations in /nl with an oscillation frequency 
~ In Kj <j)f set by the internal parameters of the theory. 

5. Conclusions 

We have studied in details by numerical and analytical methods the non-Gaussinity generated 
by a linear single field inflationary potential with superimposed oscillations. This generic form 
of the potential is an inherent feature of axion monodromy inflation. 

The model presents an example of an inflationary model embedded in string theory and can 
have a very distinct signature on observables. We show that in many cases the main feature of 
the model is a large non-Gaussianity which is detectable by future experiments such as Planck 
whereas the impact on the power spectrum is undetectably small. Thus, it appears that the 
best way to verify this type of models is through a detection of a large bispectrum with the 
very distinct shape of non-Gaussianity predicted by the model. However, very interestingly it 
also shares the usual predictions of chaotic inflation type models, i.e. a large tensor-to-scalar 
ratio and a red spectral index. 

Furthermore estimates of the trispectrum in the double-squeezed and in the counter 
collinear limit shows that non-Gaussianity from the primordial trispectrum could also be 
detectable by Planck in the present model, although a more careful study of the trispectrum 
has been left for future studies. 
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